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Abstract 

The primordial abundance of 7 Li as predicted by Big Bang Nucleosynthesis (BBN) is more than a factor 
2 larger than what has been observed in metal-poor halo stars. Herein, we analyze the possibility that this 
discrepancy originates from incorrect assumptions about the nuclear reaction cross sections relevant for BBN. 
To do this, we introduce an efficient method to calculate the changes in the 7 Li abundance produced by 
arbitrary (temperature dependent) modifications of the nuclear reaction rates. Then, considering that 7 Li is 
mainly produced from 7 Be via the electron capture process 7 Be + e~ — > 7 Li + v e , we assess the impact of the 
various channels of 7 Be destruction. Differently from previous analysis, we consider the role of unknown 
resonances by using a complete formalism which takes into account the effect of Coulomb and centrifugal 
barrier penetration and that does not rely on the use of the narrow resonance approximation. As a result of 
this, the parameter space for a nuclear physic solution of the 7 Li problem is considerably reduced. We exclude 
that resonant destruction in the channels 7 Be + t and 7 Be + 3 He can explain the 7 Li puzzle. Resonances in 
7 Be + d and 7 Be + a could potentially produce relevant effects but very favorable conditions are required. For 
the 7 Be + a channel, the possibility of a (partially) suitable resonant level in n C is studied in the framework 
of a coupled-channel model. 
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1 Introduction 



Big Bang Nucleosynthesis (BBN) is one of the solid pillars of the standard cosmological model and represents 
the earliest event in the history of the universe for which confirmable predictions can be made (see 0]] for a 
review). The theory predicts that relevant abundances of light elements, namely 2 H, 3 He, 4 He and 7 Li, were 
produced during the first minutes of the evolution of the universe. Theoretical calculations of these abundances 
are well defined and are very precise. The largest uncertainties arise from the values of cross-sections of the 
relevant nuclear reactions and are at the level of 0.2% for 4 He, 5% for 2 H and 3 He and 15% for 7 lQ 

In standard BBN, the primordial abundances depend on only one free parameter, the present baryon-to- 
photon ratio r\ = (Nb -N^)/N y , which is related to the baryon density of the universe by Q^h 2 = 3.65 • 10 7 tj. This 
quantity can be constrained with high accuracy from the observation of the anisotropics of the Cosmic Microwave 
Background (CMB). The latest WMAP-7 results suggest Q, B h 2 = 0.02249 ± 0.00056, which corresponds to 
Vcmb = 6.16 ±0.15 x 10 0. If this value is accepted, then BBN is a parameter free theory which can be used 
to test the standard cosmological model and/or the chemical evolution of the universe. 

Comparison of theoretical predictions with observational data is not straightforward. Data are subject to 
poorly known evolutionary effects and there are systematic errors. Even so, the agreement between the predicted 
primordial abundances of 2 H and 4 He and the values inferred from observations is non-trivial. However, the 
situation is much more complicated for 7 Li. Using rj = tjqmb, the predicted primordial 7 Li abundance is (6l 

(Li/H) BBN -(5.l!;j; 7 )xl0- 10 . (1) 

This is a factor ~ 3 larger than that inferred by observing the so-called 'Spite Plateau' in the 7 Li abundance of 
metal-poor halo stars, which has been given [TJ as 

(Li/H) obs a (1.7 ± 0.06 ± 0.44) x 10~ 10 . (2) 

The quoted errors take into account the dispersion of the various observational determinations. Moreover, it is 
considered that Lithium in Pop II stars can be destroyed as a consequence of mixing of the outer layers with the 
hotter interior. This process can be constrained by the absence of significant scatter in Li versus Fe in the Spite 
Plateau QQ. 

The abundance of 7 Li is a central unresolved issue in BBN [ 6 ] about which there has been recent concern (TJ 
El |3 regarding erroneous evaluation of nuclear reaction rates responsible for 7 Li production. At rj = ?7cmb > 7 Li 
is mainly produced from 7 Be via the electron capture process e~ + 7 Be — > 7 Li + v e . Thus nuclear reactions 
producing and destroying 7 Be must be considered. The leading processes, 3 He(ar, y) 7 Be and 7 Be(?i, p) 7 hi, have 
been well studied and the cross sections are known to a few percent accuracy 0. In |[T0l , it was noted that an 
increase by a factor greater than 1000 in the sub-dominant 7 Be(<f , p)2a cross section could provide the necessary 
suppression of 7 Li. This enhancement was not found in experimental data ifTTl but could have escaped detection 
if it were produced by a sufficiently narrow resonance, as suggested in [12]. Other possible resonant destruction 
channels have been considered [ 13 ] as well, such as the channels 7 Be + 3 He — > 10 C and 7 Be + 1 — > 10 B that await 
experimental verification. 

In this paper, we consider further the cosmological 7 Li problem from the nuclear physics perspective. To 
do this, in Sec. |2j we introduce an efficient method to calculate the response of the 7 Li primordial abundance 
to an arbitrary modification of the nuclear reaction rates. This approach leads to an understanding, in simple 
physical terms, of why it is so difficult to find a nuclear physics solution to the observed discrepancy. Then, in 
bees, the various possible 7 Be destruction channels are considered, including possible new resonances. 

1 After the first evaluations | 2|, theoretical uncertainties in BBN have been carefully assessed in several papers. Monte-Carlo and 
semi-analytical approaches have been used. As examples see Refs. (3]0- 
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This has been suggested recently also in Ref. Ifl3l . but here we use a more refined description of the nuclear 
processes than in previous investigations. In particular, we do not assume the narrow resonance approximation 
and we include the effect of the Coulomb and centrifugal barrier penetration in our parametrisation of resonating 
cross sections. As a result, the parameter space for a nuclear physics solution of the 7 Li puzzle is considerably 
reduced. Our conclusions are summarised in Secj6] 



2 The 7 Li response to nuclear reaction rate modifications 



It is useful to briefly review the light element production mechanism in the early universe. The abundance of a 
generic element i in the early universe, F; = n[/n B where is the baryon number density, evolves according to 
the rate equations, 

dYi 



dt 



«B 



YjYjYkio-jkvyT-J^YMo-uvfr . (3) 

j,k i 

The sums include the relevant production and destruction reactions and (cr^ v)j are the thermally averaged cross 
sections. It is known Ifl4l that a good approximation is obtained by studying the quasi-fixed point of the above 
equations, viz. 

„ Ci 



D, 



T=Tjf 



where C, and D, are the total rate of creation and destruction of the i-element, given by 

Q = n B ^ Yj Y k (o-j k v) T , 



and 



A = Yi(o-uv) T . 



(4) 



(5) 



(6) 



Ttf is the freeze-out temperature for the i-element, namely the temperature below which the rates C\ and D; 
become smaller that the Hubble expansion rate (see Ref. lTT4ll for details). 

At 77 = 7/cmb - 6 x 10~ 10 , 7 Li is mainly produced from 7 Be that undergoes at late times (i.e. long after the 
7 Be freeze-out) the electron capture process e~ + 7 Be — » 7 Li + v e , so that we have 

Y U ~ Yfte ~ yr~ , (V) 

^Be r=r B e,f 

where Cfie and D Be are the total 7 Be production and destruction rates. The dominant 7 Be production mechanism 
is through the capture reaction, 3 He(a, y) 7 Be; a reaction that has been studied in detail both experimentally lfT31l 
and theoretically |[T6ll . The cross section is known to ~ 3% uncertainty. The dominant 7 Be destruction channel 
is the process 7 Be(«, p) 7 Li; the cross section of which now is known to ~ 1% accuracy as we discuss in the next 
section. As these leading processes have been well studied, a sizeable reduction of the 7 Li predicted abundance 
occurs only if large increases of the sub-dominant 7 Be destruction channels are allowed as may be the case if 
new, so far unknown, resonances become influential. To study this possibility, we introduce a simple formalism 
to describe the response of 7 Li to a generic (temperature dependent) modification of the nuclear reaction rates. 
Motivated by Eq.(|7]), we assume that a linear relation exists between the inverse 7 Li abundance, Xu = l/Yu, 
and the total 7 Be destruction rate. This relation can be expressed in general terms a^J 

dT 



r, 

SXu-j- 



K(T) 5D Be (T), 



(8) 



2 Here and in the following, the notation Q refers to the standard value for the generic quantity Q. 



3 




where 

SX Li = =^ - 1, (9) 

and 

SD Be (T) = - 1. (10) 

D Be (T) 

The integral kernel K(T) has been evaluated numerically by considering the effects of localised (in temperature) 
increases of the reaction rate D Be (T). Our results are shown by the black solid line in the left panel of Fig. [T] 
We checked numerically that Eq. ^ describes accurately large variations of the 7 Li abundance (up to a factor 
~ 2) and, thus, it is adequate for our purposes. The kernel K(T) is peaked at ~ 50 keV; roughly corresponding 
to the 7 Be freeze-out temperature Tee/- The area under the curve is equal to - 0.7 which indicates that the total 
destruction rate of 7 Be should be increased by a factor ~ 2.5 to obtain a factor 2 reduction in the abundance of 
7 Li. 

We can use Eq. ^ to assess the sensitivity of the abundance of 7 Li with respect to a specific reaction channel. 
The total 7 Be reaction rate is given by 

DBc(T) = n B Y i Y aCn<tr a v) T , (11) 

a 

where cr a is the cross section of the reaction 7 Be +a and Y a represents the elemental abundance of the a nuclei. In 
standard BBN, the dominant contribution is provided by the 7 Be(«, /?) 7 Li reaction; accounting for about ~ 97% 
of the total 7 Be destruction rate. The standard rate D Be (T) can then be set with a few percent accuracy by, 

D Bs (T)^n B Y n {T)(a av v) T , (12) 

where Y n (T) is the neutron abundance and cr np is the cross section of 7 Be(n, p) 7 hi. Sub-dominant reaction 
channels can provide a non-negligible contribution only if there is a large increase of their assumed cross section 
values. The fractional enhancement of D Bc (T) due to a generic 7 Be + a process can be evaluated from, 

SD BeA T) = Y J^p^, (13) 
Y n (T) <o- np v>r 
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under the reasonable assumption^] that the inclusion of a new channel for 7 Be destruction does not alter the 
abundance of the a nuclei. In the right panel of Fig. [TJ we show the light element abundances, Y a (T), calculated 
assuming tj = tjcmb in standard BBN. We see that d, 3 He and 4 He have abundances larger or comparable to that 
of neutrons at the temperature T ~ 50 keV relevant for 7 Be synthesis. Thus reactions involving these nuclei could 
provide a non-negligible contribution to DB e (T) even if their cross sections are lower than that of 7 Be(w, p) 7 hi. 
This point can be expressed quantitatively by rewriting Eq. ([8]) as 



J — KJJ) t^-^-, (14) 



— K a (T) 



where _ 

Y (T) 

K a (T) = K{T)J^. (15) 

Y n (T) 

The kernels K a (T) are shown in the left panel of Fig. [I] for the cases a = n, d, t, 3 He, 4 He and can be used 
to quantify the requirements for a nuclear physics solution of the cosmic 7 Li problem. To be more consistent 
with observed data, a reduction of the 7 Li abundance by a factor 2 or more is required and that corresponds to 
SXu > 1. To obtain this, the ratios R a = (o- u v)t /(o r npv)r at temperatures T - 10 - 60keV should be, R n > 1.5 
for additional reactions in the 7 Be + n channel, Rj > 0.01 for reactions in the 7 Be + d channel, R t > 1.5 for 
reactions in the 7 Be + t channel, Ruc3 > 0.03 for reactions in the 7 Be + 3 He channel, and Ru e 4 > 4 x 10 -6 for 
reactions in the 7 Be + 4 He channel. In the next section, we explore these possibilities on the basis of general 
nuclear physics arguments. 



3 Treatment of nuclear cross-sections 

Except for spin statistical factors, the partial reaction cross section for a generic process, 7 Be+a, cannot be larger 
than cr max = (21 + \)nX 2 where / is the orbital angular momentum of the channel considered, X = l/k, and k is 
the momentum in the center of mass (CM). The relation can be rewritten as 

<r max = (2/+l)7r^ 2 = (2/+l)-^- (16) 

2/j£ 

where E is the CM energy and fi is the reduced mass of the colliding nuclei. 

For low-energy reactions involving charged nuclei and/or a non-vanishing angular momentum, the cross 
section is suppressed due quantum tunnelling through the Coulomb and/or centrifugal barrier. Modelling the 
nuclear interaction potential by a square well with a radius R, the partial cross section for the formation of a 
compound system C in the process 7 Be + a — > (C) can be expressed as, 

0"C - C"max T t . (17) 

The factor T\ represents the transmission coefficient for the channel considered. In the low-energy limit, k «; K, 
it can be calculated from, 

4k 

Ti = — v/. (18) 

K is the relative momentum of the particles inside the compound system. The functions vi are known as pene- 
tration factors iTTTl. 



3 As we see in the right panel of Fig.[TJ the abundance of Beryllium is much lower than the abundances of d, t, 3 He and 4 He. This 
implies that a tiny fraction of these elements can potentially produce a very large depletion of Beryllium nuclei. 
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For uncharged particles, the penetration factors v/ are given by 



v; 



G](R) + Fj(R) 



(19) 



where Fi(R) and Gi(R) are the regular and irregular solutions of the Schrodinger radial equation, which have 
been tabulated in Ref. [18]. Feshbach and Lax in 1948 also tabulated the associated functions v' r For the lowest 
angular momenta, 



vo = 1 



v' Q =l 



Vi = 



1 + x 2 



1 / 1 

~ 2 + \ l -x 2 



where x = kR - ^2/jE R. The transmission coefficient for reactions involving neutrons can be calculated 
exactly. One obtains 



Ti = 



AxXvi 



X 1 + (2xX + x 1 v' ] )v{ 



(20) 



where X = KR. Eq. (20 1 coincides with Eq. ( 18 1 in the low-energy limit. 

For charged nuclei, one has to rely on numerical calculations. An approximate expression for the penetration 
factors v/ can be obtained by using the WKB approximation. For the collision energy E lower that the height of 
the potential barrier, 

k,(R) 



vi 



exp 



Jr 



ki(r) dr 



where ro is the classical distance of closest approach while ki(r) is given by 



yfc Uiir) - k\ 



with 



Z a Z x e 2 1(1+1) 

Ui(r) = + . 

r 2[ir L 



(21) 



(22) 



(23) 



Eq. (21), however, does not produce accurate results when the collision energy is close to the height of the 
potential barrier and then exact expressions for the penetration factors v/ have to be used. When the Coulomb 



interaction is taken into account, Eq. ( 19 ) is still valid, but now F\(R) and G\(R) are Coulomb Functions, namely 
the regular and irregular solutions of the Schrodinger radial equation that include the Coulomb potential. Such 
functions can be evaluated numerically with standard numerical techniques. 

In the presence of an isolated resonance, the cross section for a generic process 7 Be + a— >C*— + y can 
be described by the Breit-Wigner formula, 



Tin Tr, 



cr = 



2fiE(E-E r ) 2 + r* /4' 



(24) 



where E r is the resonance energy, Y m is the width of the entrance channel, r out is the width for the exit channel 
and r tot = Tin + r out + ... is the total resonance width. The first factor in the right hand side of Eq. ( 24 1 is an 
upper limit for the cross section. Basically it coincides with cr max apart from the factor, 



to 



2J C , + 1 



(2/ fl + l)(2./ 7 + l)' 



(25) 
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that takes into account the angular momenta J a and Jj of the colliding nuclei and the angular momentum Jq> of 
the excited state in the compound nucleus. The resonance width T- m (and r out ) can be expressed as the product, 



r m = 2Pi(E,R)yf B , (26) 

where the functions Pi(E, R) describe the Coulomb and centrifugal barrier penetrations and are related to the 
penetration factors v; b)|^} 

Pi(E,R) = kRvi. (27) 

The reduced width, yr n , incorporates all the unknown properties of the nuclear interior. This has to be smaller 
that the Wigner limiting width, that is given by EUl . 



4 The 7 Be + n channel 



At T ~ Tft e f - 50keV, the 7 Be(«, p) 7 hi reaction accounts for ~ 97% of the total 7 Be destruction rate. This 
process has been very well studied. Experimental data have been obtained either from direct measurements or 
from the 7 hi(p, n) 7 Be inverse reaction (see Ref. [21] and references contained therein for details). The cross 
section for this reaction near threshold is strongly enhanced by a 2~ resonance at E x = 18.91 MeV(^] In addition, 
the data show evidence for two peaks that correspond to the (unresolved) 3 + states at 19.07 and 19.24 MeV and 
to the 3 + resonance at 21.5 MeV. The reaction rate has been determined by /^-matrix fits to the experimental data 
with uncertainties ~ 1% More refined theoretical approaches based on modern coupled-channel effective 

field theory lead essentially to the same conclusion Ell . 

The cross section of 7 Be(«, p) 7 Li reaction is extremely large. The Maxwellian-averaged cross section at 
thermal energies (cr np >7- = 3.84 x 10 4 barn E31I241 is the largest thermal cross section known in the light element 
region. At the relevant energies for 7 Be synthesis, E^ e - ?Be,f - 50keV, we have cr np (EB e ) - 9 barn; a value 
quite close to the unitarity bound value of cr max (EBe) = ^/(2/i£"Be) - 15 barn. Therein /j, is the reduced mass 
of the 7 Be + n system. The size of this cross section makes it difficult to find comparable channels for 7 Be 
destruction. 

As an alternative process, we consider the reaction 7 Be(n, a) 4 He for which no experimental data exist in the 
energy range relevant for primordial nucleosynthesis. This process normally is included in BBN calculations 
by using the old reaction rate estimate EBl to which is conservatively assigned a factor 10 uncertainty Q. The 
process 7 Be(«, or) 4 He is the second most important contribution to the 7 Be rate of destruction, accounting for 
~ 2.5% of the total. Nevertheless, due to the large uncertainty assigned, it provides the dominant contribution to 
theoretical errors in the 7 Li abundance evaluations 0. 

To obtain a factor of 2 reduction in the cosmic 7 Li abundance, the cross section of the 7 Be(«, a) 4 He reaction 
need be increased ~ 60-fold to have cr nff (2s Be ) 1.5 cr n p(^Be) ~ 15 barn; an extremely unlikely possibility. 
An upper bound on the non-resonant contribution to cr na can be obtained by considering the upper limit on the 
Maxwellian-averaged cross section {<T na )T < 0.1 mbarn that was derived [26] using thermal neutrons. Due to 
parity conservation of strong interactions, the process cannot proceed via an s-wave collision. If we assume a 
/j-wave collision, the measured value can be rescaled to r Be ,f ~ 50keV according to (tT na )T x Vr, obtaining 
the value (cr na ) TBef < 0.02 (cr np ) rBef . That result is much lower than what is required to solve the 7 Li problem. 

4 Note that in Ref. 1 19 1, the quantities Pi(E,R) themselves are referred to as the penetration factors. 

5 The entrance energy of the 7 Be + n channel with respect to the 8 Be ground state is E- m = 18.8997 MeV. Thus the listed resonances 
correspond to a collision kinetic energy equal to E r = E x - E- m = 0.01, 0.17, 0.34 and 2.6 MeV, respectively. 
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One can question this estimate because it involves extrapolation over several orders of magnitude. Irrespective 
of this, the process 7 Be(?i, or) 4 He will be suppressed at low energies with respect to 7 Be(«, p) 7 hi because of 
centrifugal barrier penetration; a quantitative estimate of which can be obtained by using the results discussed in 
the previous section. At low energy, cr nQ ,/cr np ~ T\/Tq ~ 2/uER 2 where To (T\) is the transmission coefficient 
for an s-wave (/?-wave) collision in the 7 Be + n entrance channel. By considering E = E^ e , and by taking 
R < 10 fm as a conservative upper limit for the entrance channel radius, we obtain cr na (2sBe) ^ 0.2 cr n p(^Be); 
that is also insufficient to explain the 7 Li discrepancy. 

Finally, we note that we do not expect a large resonant contribution to the 7 Be(«, a) 4 He cross section. The 
8 Be excited states relevant for 7 Be(n, p) 7 Li reaction, due to parity conservation, do not decay by of-emission. In 
summary, in view of experimental and theoretical considerations, it appears unlikely that the 7 Be + n destruction 
channel is underestimated by the large factor required to solve the 7 Li problem. 



5 Other possible 7 Be destruction channels 

In standard BBN, the 7 Be destruction channels involving charged nuclei are strongly sub-dominant. To produce 
sizeable effects on the 7 Li abundance, their efficiency has to be increased by a very large factor. This seems 
possible only if new unknown resonances are found. We discuss this possibility by using the Breit-Wigner 
formalism. 

For resonances induced by charged particle reactions at energies below the Coulomb barrier, the partial width 
of the entrance channel Ti n varies very rapidly over the resonance region. That is due to the energy dependence 
of Pi(E,R), see Eq.([26]). In general, the partial width r out of the exit channel varies more slowly since the energy 
of the emitted particle is increased by an amount equal to the Q-value of the reaction. We consider such a case 
and so neglect energy dependence of r out to have a cross section form, 

cr a = — — -z , (29) 



2»E [ ( £ _ £ r)/r 2j + [ 2Pl (E,R) + t] 2 /4 



where it has been assumed that r tot T^ + r out and £ = r out /y 2 . Then, for any chosen values of (E r , and 
for any energy E, the cross section is an increasing function of the reduced width yr . To maximise this cross 



section, we assume that the reduced width of the entrance channel is equal the Wigner limiting width (Eq. (28 )) 
that represents the maximum possible value in the approximation that the interaction potential is modelled as a 
square well of radius R. Moreover, we assume that the entrance channel is an s-wave (/ = 0) and that the factor 



<d has the maximum value allowed by angular momentum conservation by setting Jc* - J a + Jse in Eq.(25i. 
Under these assumptions, the cross section of the resonant process is 

a a = — — -j — . (30) 

^E YE-E r )ly%(R)\ +[2P (E,R) + t] 2 /4 

It is a function of the two resonance parameters, (E r , and of the entrance channel radius R. 



We have applied Eq. ( 30 1 to a generic 7 Be destruction channel involving charged nuclei and then determined 



the effect on the 7 Li abundance by using Eq. ( 14 1. The thermally averaged cross section (cr a v)r has been evalu- 
ated numerically without using the narrow resonance approximation. Our results are shown in Fig.[2]as a function 
of the resonance parameters (E r ,g). The 'coloured' lines represent the iso-contours for the 7 Li abundance, 

5Y U = 1-5^. (31) 
Y U 
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Figure 2: 77ie coloured lines show the fractional reduction of the primordial 7 Li abundances that can be achieved by a 
resonance in the 7 Be + a reaction. The various panels correspond to a — d, t, 3 He and 4 He, respectively, starting from the 
upper left corner. The black dashed lines correspond to the condition r lot (E r ,R) — 50keV, which is the limit for narrow 
resonance. 
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The various panels correspond to the processes 7 Be + a with a - d, t, 3 He and 4 He respectively, starting from 
the upper left corner and proceeding clockwisd^ In our calculations, we assumed the entrance channel radius to 
be R = 10 fm. That is quite a large value considering the radii of the involved nuclei but it has been chosen to 
provide a conservative upper estimate of the resonance effects. 

It is important to note that there is a maximum achievable reduction (6Yu)ma* of the 7 Li abundance for each 
reaction channel considered; a point not discussed in previous analyses lfT3l [T21 . Those analyses included the 
effects of Coulomb barrier penetration a posteriori and used the narrow resonance approximation to calculate 
(o~ a v)T- The use of the narrow resonance approximation outside its regime of validity can lead to severe over- 
estimation of the resonance effects. Indeed, using that assumption one has (o~ a v)T ~ T e ff (yur) _3/2 exp(-E r /T) 
where the effective resonance width is defined as r e ff = (ri n r out )/r tot . This expression does not predict any 
upper limit for (cr a v)T as a function of T e ff. But a limit should exist. That can be understood by considering 
o~ a < {nod) / (2fiE) for any possible choice of the resonance parameters as per Eq.( 24 1. In fact, the correct scaling 
for broad resonances is given by {o- a v)j ~ (r e ff/r tot ) yT 1 ^ 1 T~ 3 ^ 2 where the factor r e ff/r tot = (ri n r out )/T 2 ot , 
cannot be larger than one. In Fig. [2] we display with the 'black' dashed line the result found using the condition 
r tot (E r , R) = Ti n (E r ,R)+r out = 50keV. The turning point of the line corresponds to the situation Ti n (E r , R) ~ T out 
and it is localised in the region where £ ~ 2 Po(E r , R). Results with a narrow resonance approximation can only 
be compared with those shown in the lower left corner of the plots where Yt ot (E r ,R) «: E r , and the 7 Li reduction 
typically is negligible. 

The results that we have obtained for each specific reaction channel are now outlined sequentially: 



7 Be + d: With this initial channel, the maximum achievable effect is a ~ 40% reduction of primordial 7 Li 
abundance. This reduction could substantially alleviate the discrepancy between theoretical predictions and 
observational data which differ by a factor ~ 3 in the standard scenario. The maximal effect is obtained for 
a resonance energy E T ~ 150 keV with a total width r tot (E r ,R) ~ 45keV and partial widths approximately 
equal to r out ~ 35keV and Ti n (E r ,R) ~ lOkeV. The dependence of the maximum reduction (5Fu)max on the 



assumed entrance channel radius R is shown in the left panel of Fig. 2 in which (6Yu)ma\ increases with R. 
That is expected given that the partial width of the entrance channel is determined primarily by Coulomb barrier 
penetration as in Eq. (26 1. So quite large values for R are needed to solve the cosmic 7 Li problem. Even if these 
are much larger than the sum of the radii of the involved nuclei, 2.65 fm and 2.14 fm for 7 Be and the deuteron 
respectively, they cannot be excluded in view of the large uncertainties in the approach. Our results basically 
coincides with those found [12J using a different approach. Note that there is a (poorly known) excited state 
in 9 B at 16.71 MeV excitation. It lies just 220 keV above the 7 Be + d threshold and it decays by gamma and 
particle (proton or 3 He) emission. The process 7 Be + d — > 6 Li + 3 He could enhance 6 Li production in the early 
universe. However, the effect on the final 6 Li abundance is definetely too small to explain the observations of 

Ref.rg. 



7 Be + t: With this initial channel, the maximum achievable effect is a ~ 0.2% reduction of primordial 7 Li 
abundance. The existence of a resonance in this channel cannot solve the cosmic 7 Li problem, since to produce 
a significant 7 Li reduction, the 7 Be destruction rate due to the 7 Be + 1 reaction should be comparable to that form 
7 Be + n processes. Clearly that is impossible because: 

1 . neutrons are more abundant than tritons at the relevant temperature T# e ~ 50 keV as is seen in Fig. [T] 

6 We do not consider the 7 Be + p entrance channel since this is known to be subdominant, see e.g. 1 13], and it is well studied at low 
energies due to its importance for solar neutrino production. 

7 The resonance parameters that maximise the 7 Li suppression are slightly dependent on the assumed radius. 
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2. the cross section of 7 Be(«, p) 7 Li is close to the unitarity bound while the 7 Be+? collisions are suppressed 
by Coulomb repulsion. 

Conclusions reached by others [13] differ from ours. Theirs are artifacts from using of the narrow resonance 
approximation outside of its regime of application. 

7 Be + 3 He: With this channel, the maximum achievable effect is a ~ 10~ 4 reduction in the abundance of 
primordial 7 Li. Again then, a resonance in this channel cannot solve the cosmic 7 Li problem. The small effect 
is due to strong Coulomb repulsion suppressing the partial width of the entrance channel (Ti B (E r , /?)). From 
our calculations, we obtain Ti n (E r ,R) ~ 0.2 meV and T{ n (E r ,R) ~ 7eV for E r - 100 keV and E r - 200 keV 
respectively; values that are much smaller than required to obtain non negligible effects. 

7 Be + a: With this channel, the maximum achievable effect is a ~ 55% reduction of primordial 7 Li abun- 
dance. This is obtained for a resonance with a relatively large centroid energy E r ~ 360 keV, with a total width 
T tot (E r ,R) ~ 21 keV and with partial widths T out ~ 19keV and T m (E r ,R) ~ 1.5 keV. The strong suppression of 
the cross section due to Coulomb repulsion in this case is compensated by the fact that the a nuclei are ~ 10 6 
times more abundant than neutrons when the temperature of the universe falls below ~ 70keV. But this is limited 
to specific regions of the parameter space, onl>[^] The dependence of the maximal achievable reduction (<5lLi)max 
from the entrance channel radius R is explored in the right panel of Fig. [3] We see that R > 9 fm is required to 
obtain (SYu)max ^ 0-5 and, thus, to have a substantial impact on the 7 Li puzzle. 

The possibility of a missing resonance in the last channel is particularly intriguing theoretically. We have 
calculated the spectrum of U C using a coupled-channel model for the n- l0 C system, with coupling involving the 
excited quadrupole state of 10 C. A multi-channel algebraic scattering (MCAS) was used with which account is 
made of constraints imposed by the Pauli principle on single -particle dynamics besides coupling interactions to 
the collective excitations of the 10 C states E71l . Bound and resonant low-energy spectra of light nuclei have been 
analysed systematically with this approach |[28l|29l and in particular for carbon isotopes f30l I3T1 1321 . In Table[T] 
a comparison is given between the calculated spectrum with the observed levels of n C [33]. Clearly there is a 
one to one correspondence except for a \ state predicted at 6.885 MeV. That excitation energy lies relatively 
close to entrance of 7 Be + 4 He channel which is 7.543 MeV above the U C ground state and would require a 
cf-wave collision (or a coupled-channel transition to the 7 Be first excited state) to ensure angular momentum 
and parity conservation. It is worthwhile therefore to investigate this system further. To find if there is a (so far 
missing) resonance with about the right energy that could significantly increase cross sections. 

6 Conclusions 

We have investigated the possibility that the cosmic 7 Li problem originates from incorrect assumptions about 
the nuclear reaction cross sections relevant for BBN. To do so, we introduced an efficient method to calculate 
the changes in the 7 Li abundance produced by an arbitrary (temperature dependent) modification of the nuclear 
reaction rates. Then, taking into account that 7 Li is mainly produced through 7 Be, we used this method to 
assess whether it is possible to increase the total 7 Be destruction rate to the level required to solve (or alleviate) 

8 In Fig. 9 of Ref. 1131 . it was shown that a narrow resonance at an energy E < 100 keV and effective width r cff ~ lOmeV can 
substantially reduce the 7 Li abundance. According to our analysis, however, this region of parameter space is unphysical since, if we 
take into account correctly the effects of Coulomb repulsion, we obtain Ti n (E r ,K) < 1 meV for E r < 100 keV. 
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J p Nuclear data MCAS levels 



3/2" 
1/2- 
5/2" 
3/2" 
1/2+ 
7/2- 
5/2+ 
1/2" 
3/2+ 
3/2" 
5/2" 
7/2+ 
5/2+ 
5/2+ 
3/2" 
5/2" 
7/2" 



0.00 
2.000 
4.3188 
4.8042 
6.3392 
6.4782 
6.9048 



7.4997 

8.1045 

8.420 

8.655 

8.699 

9.20 

9.65 

9.78 

9.97 
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0.00 

2.915 

3.225 

3.303 

8.373 

5.768 

7.781 

6.885 

11.059 

7.332 

9.689 

10.343 

10.698 

11.868 

11.253 

12.802 

9.022 



Table 1: Spectra of C. The data are taken from Ref. \33\l while calculated values have been obtained with the MCAS 
formalism [27]. Potential parameters defining the MCAS coupled-channel interactions have not been sought for the optimal 
reproduction of levels, but only to check for possible missing resonances. 

the cosmic 7 Li puzzle. Given present experimental and theoretical constraints, it is unlikely that the 7 Be + n 
destruction rate is underestimated by the ~ 2.5 factor required to solve the cosmic 7 Li problem. 

On the basis of very general nuclear physics considerations, we have shown that the only destruction channels 
that could have a relevant impact on the 7 Li abundance are 7 Be + d and 7 Be + a. Our results suggest that it is 
unrealistic to consider new resonances in 7 Be + ? and 7 Be + 3 He channels to solve the 7 Li problem. However, with 
the other channels, new resonances must exist at specific energies and with suitable resonance widths. Postulating 
a resonance in the 7 Be+<i reaction at an energy E r ~ 150 keV, with a total width T tot (E r , R) ~ 45 keV, and partial 
widths r out ~ 35 keV and ri n (E r , R) ~ 10 keV gave a ~ 40% reduction in the 7 Li abundance. A larger suppression 
of ~ 55% was obtained by assuming a resonance in the 7 Be + 4 He channel with energy E r ~ 360 keV, with a 
total width T tot (E r ,R) ~ 21 keV and partial widths r ou t ~ 19 keV and T{ n (E r ,R) ~ 1.5 keV. These results are the 
maximal achievable reductions of the 7 Li abundance, since they were obtained assuming that the resonance width 
of the entrance channel has the largest value allowed in the presence of Coulomb repulsion, and by scanning the 
space of the other resonance parameters. Also, we considered a relatively large value for the entrance channel 
radii,/? - 10 fm. 

In summary, a nuclear physics solution of the 7 Li puzzle requires very favourable conditions. Due to the 
importance of the problem, and before considering more exotic explanations, further experimental and theoretical 
analysis of the 7 Be + d and 7 Be + a reaction channels, specifically, are necessary. Other 7 Be destruction channels 
do not appear influential. The 7 Be + a reaction is of particular interest because the comparison of the observed 
and calculated spectra for n C have one-to-one correspondence except for a predicted I state that lies relatively 
close to the energy range of interest. 
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R (fin) R (fm) 

Figure 3: The maximal fractional reduction (5FLi)max of the Lithium abundance that can be achieved by a resonance in 
7 Be + d and 7 Be + a channels as a function of the assumed entrance channel radius R. 

References 

[1] K. Nakamura et al. [Particle Data Group], J. Phys. G 37, 075021 (2010). 

[2] M. S. Smith, L. H. Kawano and R. A. Malaney, Astrophys. J. Suppl. 85, 219 (1993); 
L. M. Krauss and P. Romanelli, Astrophys. J. 358, 47 (1990). 

[3] P. D. Serpico, S. Esposito, F. Iocco, G. Mangano, G. Miele and O. Pisanti, JCAP 0412, 010 (2004). 

[4] G. Fiorentini, E. Lisi, S. Sarkar and F. L. Villante, Phys. Rev. D 58, 063506 (1998); 
E. Lisi, S. Sarkar and F. L. Villante, Phys. Rev. D 59, 123520 (1999); 
R. H. Cyburt, Phys. Rev. D 70, 023505 (2004); 



P. Descouvemont, A. Adahchour, C. Angulo, A. Coc and E. Vangioni-Flam, arXiv:astro-ph/0407101 

[5] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 192, 18 (2011). 

[6] R. H. Cyburt, B. D. Fields and K. A. Olive, JCAP 0811, 012 (2008). 

[7] S. G. Ryan et al., Astrophys. J. 530, L57 (2000); 
L. Sbordone et al., A&A 522, 26 (2010); 
P. Bonifacio et al., Astron. Astrophys. 390, 91 (2002); 
J. Melendez and I. Ramirez, Astrophys. J. Lett. 615, 33 (2004). 

[8] S. G. Ryan et al., Astrophys. J. 523, 654 (1999). 

[9] M. Asplund, D. L. Lambert, P. E. Nissen, F. Primas and V. V. Smith, Astrophys. J. 644, 229 (2006) 

[10] A. Coc, E. Vangioni-Flam, P. Descouvemont, A. Adahchour, C. Angulo, Astrophys. J. 600, 544 (2004). 

[11] C. Angulo et al, Astrophys. J. 630, L105 (2005). 

[12] R. H. Cyburt and M. Pospelov, |arXiv:0906.4373l [astro-ph.CO] . 

[13] N. Chakraborty, B. D. Fields and K. A. Olive, Phys. Rev. D 83, 063006 (201 1). 

[14] R. Esmailzadeh, G. D. Starkman, S. Dimopoulos Astrophys. J. 378, 504 (1991). 



13 



[15] F. Confortola et al. [LUNA Collaboration], Phys. Rev. C 75, 065803 (2007) [Phys. Rev. C 75, 069903 
(2007)]. 

[16] L. Canton and L. G. Levchuk, Nucl. Phys. A 808 192 (2008). 

[17] J. M. Blatt, V. F. Weisskopf, "Theoretical Nuclear Physics", Dover Publications (October 18, 2010). 

[18] H. A. Antosiewicz (Ch. 10), M. Abramowitz (Ch.14) "Handbook of Mathematical Functions", eds. M. 
Abramowitz and I. Stegun; (Nat. Bur. Stds., New York, 1964) 

[19] C. Iliadis, "Nuclear Physics of Stars", Wiley- VCH; 1 edition (April 24, 2007). 

[20] T. Teichmann and E.P Wigner, Phys. Rev. 87, 123 (1952). 

[21] A. Adahcour and P. Descouvemont J. Phys. G. 29, 395 (2003). 

[22] V. Lensky and M. C. Birse, Eur. Phys. J. A 47, 142 (2011). 

[23] F. Ajzenberg-Selove, Nucl. Phys. A 490, 1 (1988). 

[24] P. E. Koehler et al. Phys. Rev. C 37, 917 (1988). 

[25] R. V. Wagoner et al, Astrophys. J. 148, 3 (1967); 
R. V. Wagoner, Astrophys. J. Suppl. 18, 247 (1969); 
R. V. Wagoner, Astrophys. J. 179, 343 (1973). 

[26] P. Bassi et al., Nuovo Cimento 28, 1049 (1963). 

[27] K. Amos, L. Canton, G. Pisent, J. P. Svenne and D. van der Knijff, Nucl. Phys. A 728, 65 (2003). 

[28] L. Canton, G. Pisent, J. P. Svenne, D. van der Knijff, K. Amos and S. Karataglidis, Phys. Rev. Lett. 94, 
122503 (2005). 

[29] L. Canton, G. Pisent, J. P. Svenne, K. Amos and S. Karataglidis, Phys. Rev. Lett. 96, 072502 (2006). 

[30] J. P. Svenne, K. Amos, S. Karataglidis, D. van der Knijff, L. Canton and G. Pisent, Phys. Rev. C 73, 027601 
(2006). 

[31] S. Karataglidis, K. Amos, P. Fraser, L. Canton and J. P. Svenne, Nucl. Phys. A 813, 235 (2008). 

[32] K. Amos, L. Canton, P. Fraser, S. Karataglidis, J. Svenne and D. van der Knijff, "Linking the exotic structure 
of 17 C to its unbound mirror 17 Na," Nucl. Phys. A (2012) DOI:10.1016/j.nuclphysa.2012.01.022 

[33] F. Ajzenberg-Selove Nucl. Phys. A 506, 1 (1990). 



14 



